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We report high resolution polarized infrared studies of the quarter-filled spin ladder compound 
o/-NaV205 as a function of temperature (5 K< T < 300 K). Numerous new modes were detected 
below the temperature T c = 34 K of the phase transition into a charge ordered nonmagnetic state 
Q_) ■ acconrpanied by a lattice dimerization. We analyse the Brillouin zone (BZ) folding due to lattice 

dimerization at T c and show that some peculiarities of the low-temperature vibrational spectrum 
come from quadruplets folded from the BZ point ( = , 5, ?)■ We discuss an earlier interpretation of the 
70, 107, and 133 cm -1 modes as magnetic bound states and propose the alternative interpretation 
as folded phonon modes strongly interacting with charge and spin excitations. 
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PACS numbers: 78.30.-j, 63.20.-e, 63.20.Ls 

o 

O ; I. INTRODUCTION 

The magnetic insulator a'-NaV^Os has received a considerable attention since 1996, when Isobe and Ueda reported 
£> ■ a spin-Peierls-like transition at T c = 35 K H. The structure of this vanadate consists of two-leg ladders formed by 
0^ ' corner-sharing distorted V2O5 pyramids running along the orthorhombic b-axis, with their rungs along the a-axis 
(see Fig. 1). Neighboring ladders are linked via common edges of the pyramids to form the a6-layers. Na atoms lie 



between the layers. At room temperature, the structure is centrosymmetric (space group Pramn-D)^) with only one 
crystallographic position for vanadium having a formal valence +4.5 @-0]. As a result, one electron is shared by two 



o 



V atoms and is actually distributed over a V-O-V molecular orbital at the ladder rung. A quarter-filled ladder along 
the 6-direction is formed, which explains one-dimensional magnetic properties of a'-NaV^Os above T c P||. 

Below Tc, a steep isotropic decrease of the magnetic susceptibility corresponding to a singlet-triplet gap of 
A=85 K HII and a lattice dimerization described by the propagation wave vector k = (555) Q were recently found 
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to be accompanied by a charge ordering of a zigzag type |W,11-13 The nature of the state below T c is presently under 
intense investigation. The recent synchrotron X-ray diffraction studies of the low-temperature (LT) structure p, [bt|Jl5[ 
resulted in the acentric space group Fmml-C^ with three essentially different positions for the vanadium atoms ar- 
ranged into a sequence of zigzag modulated ladders separated by nonmodulated ones (see Fig. 1). The corresponding 
three vanadium valence states, +4, +5, and +4.5, have been suggested from the bond-valence analysis in Ref. p6|. 
This finding contradicts the NMR observation of only two inequivalent electronic states of vanadium below T c [|7|. 
More recent bond-valence calculations suggest, however, only two vanadium valences despite the three crystallograph- 
ically inequivalent V sites p5[ . Nevertheless, several observations, in particular, the properties of X-ray anomalous 
scattering [ p.3| , can be better explained if one assumes the modulation of all the vanadium ladders, and therefore dif- 
ferent subgroups of Fmm2 were considered J17],|l8||. In their Raman scattering study of folded modes, the authors of 
Ref. Q claimed the centrosymmetric C2/c group for the LT-phase of a'-NaV^Os, from purely spectroscopic reasons. 
Folded modes of the dimerized lattice of ev'-NaV^Os have been observed first in the far infrared (FIR) transmis- 
sion pfl|, then Raman scattering 0,B^,E2] and FIR reflection P3Lp4j measurements. Only two out of more than thirty 



observed [E0pq,E3,E|l new low-temperature FIR modes, namely, the modes at 718 cm -1 (E || a, b) p3,E6 24 and 
960 cm~ (E || a) |p4| have been studied in detail. Three lowest frequency Raman lines (67, 107, 133 cm - " that 
emerge below T c have been attributed to magnetic singlet bound states ^,^^8). The recent work [g{| called this 
interpretation into question, however. Thus, both the low-temperature structure of ev'-NaV^Os and the nature of 
the low-energy excitations still remain controversial. The necessary steps towards solving these problems include (i) 



an acquisition of precise high resolution experimental data on the properties of new low-temperature infrared and 
Raman modes in all the principal polarizations and (ii) a theoretical analysis of the properties of modes folded into 
the Brillouin zone (BZ) center of the LT-phase from different BZ-points of the paraphase. 

The present work is devoted to the detailed high-resolution FIR transmission study of the frequencies, widths, and 
shapes of additional low-temperature modes of a'-NaV20s for T < T c in E || a, E || b, and E |j c polarizations. To 
understand, at least qualitatively, polarization properties, splittings, and temperature behavior of the folded modes, we 
perform lattice dynamical calculations of the phonon dispersion curves throughout the Brillouin zone of the paraphase 
and of frequency changes due to lattice distortion and charge ordering below T c . Possible manifestation of magnetic 
excitations and of the interaction between lattice, charge, and spin excitations in the FIR spectra are discussed. 

II. EXPERIMENT 

Samples of stoichiometric a'-NaV20s used in this study were grown by a melt growth method using NaV03 as a 
flux P0|| - The process resulted in single crystals with dimensions from 1 x 5 x 0.5mm to 3 x 17.3 x 1.6mm along 
the a-, &-, and c-axes, respectively. Crystals from different batches were used. For transmission measurements in 
the k || c, E || a or E || b configurations, we prepared four thin samples cleaved perpendicular to the c-axis. Their 
thicknesses were 111, 45, 14, and 6/xm. For the measurements in the k || a, E || c or E || b geometries, two samples 
were used. One of them was 1.3 mm thick along the a-axes (sample 2). The other one was covered with epoxy and 
then polished to the thickness of 150 (im (sample 1). The samples were checked by X-ray diffraction, magnetization, 
and ESR measurements. They exhibited a sharp transition at about 34 K. 

The transmittance spectra in a polarized light incident normally to the sample surface were measured in the spectral 
range 25-2000 cm -1 at a resolution 0.05-1.0 cm" 1 using a BOMEM DA3.002 Fourier transform spectrometer equipped 
with a helium- vapor optical cryostat. Samples were mounted to a specially constructed cryostat inset that partially 
compensated thermal variations of the cryostat length. The temperature of the sample varied from 300 to 6 K and was 
stabilized with the precision of ±0.1 K. At each temperature, a reference spectrum was taken. For the case of sample 1 
measured in E || c polarization, transmittance is given relative to the layer of epoxy having the same thickness as 
the sample. Precision of the experimental line position was 0.05-0.5 cm -1 , depending on a particular line. Precise 
absolute wave number scale is an intrinsic property of Fourier transform spectroscopy. 

III. RESULTS 

Fig. H shows FIR polarized transmittance in the spectral range from 55 to 350 cm" 1 for a'-NaV^Os at temperatures 
above and below T c . At T ~ 40 K> T c in the displayed frequency range we found two phonon lines (ujto = 178 
and 225 cm" 1 ) in E || b polarization (for the discussion of the 215 cm -1 line see Ref. ||); three strongly asymmetric 
phonon peaks at 90, 140, and 254 cm -1 visible on a strong absorption background in E || a polarization; two strong 
phonon lines with ujto = 163 and 180 cm" 1 , lolo — 169 and 215 cm" 1 (compare with u>to — 162 and 179 cm" 1 , 
ujlo = 165 and 212 cm" 1 as determined from FIR reflectance spectra at room temperature H) and an absorption peak 
at 282 cm" 1 in E || c polarization. This peak can be assigned to a missing B\ u phonon placed at 298 cm" 1 by lattice 
dynamics calculations for the high-temperature (HT)-phase || and not observed in FIR reflectance spectra |y.|2(J . 

A broad absorption band centered at about 320 _1 is present in E || a spectrum (see also p0| , |3l| , |32"| | ) . One can 
clearly see it in Fig. 2a where the transmittance of a 14 /im thick sample is almost zero near this frequency No 
broad-band absorption of a comparable intensity exist it E || b (Fig. 2b) and E || c (Fig. 2c) polarizations (note that 
Figs. 2a and 2b refer to the 14 /im thick sample, wile Fig. 2c represents the spectra of a 150 /xm thick sample). The 
a-polarized FIR absorption continuum extends to about 1500 cm" 1 and consists of two bands centered at 300 and 
1100 cm - 1 ||J||5§. 

Drastic changes occur in the spectra of NaV205 at the transition temperature. First, the long- wavelength part 
of the broad-band absorption in E || a polarization decreases markedly, the sample becomes more transparent (see 
Fig. 2a). The changes of transmittance near the 178 cm" 1 phonon line in E || b polarization are caused by a 
sharp narrowing of this phonon below T c . A special study is necessary to clarify the question, whether the changes 
of background transmittance observed in E || c polarization can be explained by a similar reason, namely, by the 
changes of regular phonons 162, 179, and 282 cm -1 below the transition temperature, or there is a weak broad- 
band absorption that diminishes below T c . The spectra modeling (see below) has shown that, indeed, a weak broad 
absorption band that bleaches below T c is present in E || c polarization. The measured absorption coefficient at the 
frequency 115 cm" 1 is k(lo = 115 cm" 1 ) = 20 cm" 1 . To compare, for the E || a polarized broad-band absorption 
k(uj = 115 cm" 1 ) = 400 cm" 1 . 



Second, the shape of asymmetric peaks at 90 and 140 cm -1 in E || a spectrum changes to symmetric. Third, 
numerous new spectral features appear below T c in E || a, E || b, and E || c spectra. The first two phenomena were 
already discussed in our earlier publication [EOj. We have attributed a broad band in E || a polarization to two-spinon 
absorption, the asymmetric line shape being due to Fano-type resonance between lattice vibrations and a two-spinon 
continuum. The opening of a gap in the magnetic excitation spectrum below T c leads to a vanishing of a low-frequency 
continuum and to the appropriate transformation of the line shapes. Subsequent structure redetermination for the 
HT phase of a'-NaV205 ||~|| does not, however, allow to ascribe E || a polarized spectral continuum to a pure two- 
spinon absorption. Symmetry considerations show that for both HT- and LT-structure of a'-NaV20s the two-spinon 
(two-magnon) absorption is allowed only in E || c polarization [^3|,|4). Probably, a weak FIR absorption background 
observed in this polarization results just from the two-spinon absorption. In that case, the observed diminishing 
of this absorption in the low- frequency region of E || c spectra would correspond to the opening of a gap in the 
magnetic excitation spectrum below T c and also one would expect to observe the magnetic bound states, if any, just 
in E |j c polarization. As for the E || a polarized low-frequency continuum, it could be formed by two spinons and 
one low-energy charge excitation, as has been proposed recently by Khomskii et al. in their paper on the spin-isospin 
model of a'-NaV 2 5 ||. 

In what follows, we concentrate our attention on the spectral lines that appear below T c . 



A. New spectral lines 

Numerous new spectral lines appear in transmittance spectra below T c ~ 34 K and grow in intensity upon cooling 
the sample. For the spectral region 50-350 cm -1 they are marked by arrows in Fig. 2. Figs. 3-5 show some examples 
of the E |j a, E || b, and E || c absorbance difference spectra a(T) — a(37K) for several temperatures. Here, the 
absorbance is defined by a — — log 10 Tr/d, where Tr is the transmittance of a d cm thick sample. 

The temperature behavior of the majority of new lines is characterized by an appreciable broadening and small 
(0.2-2 cm' 1 ) but clearly detectable shifts near T c (see Figs. ||, 0). It should be mentioned that the vibrational lines 
that exist both below and above T c also broaden and shift when approaching T c from below (see insets of Fig.Jfl). 
Positions, widths, and oscillator strengths of the observed spectral lines that appear below T c are listed in Table ^ 

These parameters were obtained from the following fit ting procedure. The experimental transmittance spectrum 
at normal incidence was approximated by the expression |3q| 



Tr = 



4y^ 



(x/i+l) 2 e- 



-if3d _ 



(^-1) 2 < 



il-Sd 



(1) 



f3 = 2TTLUy/e 

in the case of clearly visible interference pattern, or by the expression 
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in the absence of interference within the sample. Here 
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is the absorption coefficient. The dielectric function e(co) was represented as a sum of contributions from N nonin- 
teracting oscillators 
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where e^, is the high-frequency dielectric constant, u)j,^j, and fj are the frequency, the damping constant, and the 
oscillator strength of the j-th oscillator. First of all, the calculated spectrum was fitted to the experimental one at 



T=37K, just above the temperature of the phase transition, starting from the room temperature values determined 
earlier in our work J5|. The fitting to the low-temperature (T ~10K) spectrum was performed in two steps. First, 
changes of the phonons that exist both in HT- and LT-phases were taken into account and, then, new oscillators were 
introduced. 

It should be noted that we could not study several spectral regions because of a too strong absorption of either the 
sample or mylar windows of our cryostat. In particular, the lines at 718 cm -1 are only partially seen, being too close 
to a strong absorption band of mylar. Results of the E || c, k || a measurements refer to the 150 /xm thick sample. 
We were able to detect the lowest frequency lines also in the 1.3 mm thick sample. Their positions are indicated in 
brackets in Table fil 

Fig. 8 shows the temperature dependence of the normalised oscillator strength of the modes 101.4 and 111.7 cm -1 
(E || b polarization); 126.7 cm' 1 (E || a); 70, 106.9, 124.7, 133.1, and 140 cm' 1 (E || c). The data on the 101.4 and 
111.7 cm" 1 modes summarize a very detailed study of two samples with the thicknesses d\ = 14 /im and d 2 = 111 /xm. 

The intensity of the 101.4, 126.7, and 140 cm -1 modes follows the same temperature dependence as the intensity 
of new X-ray reflections corresponding to the doubling (a, b) and quadrupling (c) of the lattice constants |l(J (also 
shown in Fig. 8). This dependence is characterized by a steep rise just below T c and a saturation below ~ 22 K. 
The majority of the observed new lines behave like that, while several ones (in particular, the lines 91.2, 145.0, 145.7, 
157.2, and 451.1 cm" 1 in E || a, 111.7 cm" 1 in E || b, 70, 106.9, 126.7, and 133.1 cm" 1 in E || c) demonstrate a more 
slow growth of intensity below T c and a saturation at lower temperatures. We could not measure reliably the intensity 
dependence of the very weak feature at about 133 cm" 1 in E || c polarization but it seems to grow in the slowest 
way among all the studied modes and to follow the same temperature dependence as the spin-gap mode observed in 
optical ESR measurements |37J (also shown in Fig. 8). 

B. Factor-group analysis (Fmm2). Brillouin zone folding for different proposed low-temperature structures. 

According to the recent X-ray diffraction work |q jl4|jl5|l , the structure of the low-temperature phase of c/-NaV2C>5 
is face-centered, with Fmm2 space group on the 2a x 2b x Ac supercell. The primitive cell corresponding to the 
translation vectors (a, 0,2c), (0,6,2c), and (a, b, 0) contains 64 atoms and, hence, there are 189 optical vibrational 
modes. Table O lists the Wyckoff positions of atoms in the Fmm2 supercell H and in the Pmmn basic structure [0. 
The corresponding local symmetries are also indicated. These positions yield the following irreducible representations 
of the vibrational modes at the Brillouin zone center of the Fmm2 LT-structurc |B8| : 
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3A X + 3A 2 + iB 1 + W 2 
A x + B x + B 2 
2A X +A 2 + 2B X + B 2 
2A X +A 2 +B X + 2B 2 
A x +A 2 + 2Bx +2B 2 . 



Multiplying the representations given above by the number of different positions of the appropriate symmetry, sum- 
ming them, and subtracting acoustic modes (Ai + B\ + B 2 ) , we obtain the following optical vibrational modes: 

r N. b v 2 o 5 ( fmm2 ) = 55Ai(aa,66,cc;£ || c) + 36 A 2 {ab) + 55Bi {ac; E \\ a) + 43B 2 (bc; E \\ b). 

153 modes (A\, B±, B 2 ) are both Raman and IR active, while 36^2 modes are allowed only in Raman scattering. 

More information on the transition-induced vibrational modes, in particular, on the distribution of their frequencies, 
can be obtained from the analysis of the BZ folding at the structural phase transition. We have constructed the 
Brillouin zones for the Pmmn HT-structure (it is shown in Fig. 9), for the Fmm2 LT-structure, and also for the other 
proposed LT-structures and analysed what points of the BZ of the paraphase fold into the new zone center and thus 
become optically active. Results of this analysis are represented in Table H. 

In the case of Fmm2 low-temperature structure, the two Z-points (0, 0, ±^) of the paraphase and the eight Q-points 
of the type (±, \, \) deliver 48 singlets 16^i + 6A 2 + 165i + 10^ (compatible with (8Zi + 8Z 2 ) + {3Z 3 + 3Z 4 ) + 
(8Z5 + 8Zq) + (5Zj + 5Z$) irreducible representations of the symmetry group of the Z-point) and, respectively, 24 
quadruplets 12(2Ai + 2A 2 ) + 12(2i?i + 2B2) (compatible with Qi and Q 2 doublets of the Q-point) to the BZ center 
of the LT phase. They join 15Ai + 6A 2 + 15i?i + 9B 2 singlet vibrational modes that originate from, correspondingly, 
{7B lu + 8A g ) + (3A U + 3B lg ) + (7B 3u + 8B 2g ) + {AB 2u + 5B 3g ) vibrations of the BZ center (r-point) of the paraphase. 
Because of a degeneracy of the Q-point folded modes, the number of additional frequencies observed below T c should 
be not so high as it followed from a simple factor-group analysis (namely, 72 instead of 189 — 45 = 144). 



In the case of the face-centered structures Ccc2, C2/c, or C2, there are 128 atoms in the primitive cell and, thus, 
381 optical modes. In addition to the just mentioned BZ points Z and Q, the A, S, and R points also fold into 
the zone center delivering 96 additional doublets. A lattice distortion in the LT-phase lifts a degeneracy of folded 
quadruplets or quadruplets and doublets and causes their splitting. We shall consider the polarization properties of 
the split components in the section IVA. 

C. Lattice dynamics calculations 

We consider the lattice dynamics of a'-NaX^Os in the framework of the rigid-ion model. The model incorpo- 
rates long-range Coulomb interactions and short-range interactions described by Born-Mayer potentials. The model 
parameters were obtained in Ref. [gj by fitting to the HT-phase zone-center frequencies measured by optical tech- 
niques. In the present work, we use them to calculate the phonon dispersion curves throughout the Brillouin 
zone of the paraphase. As an example, Fig. 10 displays the calculated dispersion curves along the closed path 
T(OOO) - A(00|) - Z(00i) - R{\\ §) - Q{1 \\) - S(iiO) - r(000) in the BZ of the high-temperature Pmmn structure 
for the two highest-frequency Davydov doublets B% u — T5i g and B\ u — A g originating from the stretching vanadium- 
apical oxygen vibration. We do not show here the rest of 48 phonon branches. Most of them exhibit pronounced 
slopes and complicated interactions. Many (but not all) of the dispersion curves along the T—Z direction are rather 
flat, reflecting more weak forces between layers in comparison with those inside a layer. However, only 12 out of 45 
optical branches show a dispersion less than 6 cm -1 (6 cm" 1 is a maximum line-width among the observed folded 
modes), and so the phonon spectrum of the LT-phase cannot be described in the framework of one layer, as has been 
proposed in Ref. |19| . 

Our lattice dynamical model neglects electronic polarization and does not take into account three-body forces 
that strongly affect the frequencies of bending vibrations. Because of that, there are several large discrepancies 
(50-80 cm -1 ) between the calculated and observed zone-center phonon frequencies of the high-temperature Pmmn 
structure [|5j and one cannot expect a good quantitative description of phonon frequencies throughout the Brillouin 
zone. (It should be mentioned, however, that the sound velocities determined from the slope of the acoustic branches 
along the ladders (the 6-axis) at k = 0, v = ($£)fe=o, v t = 7.0 ■ 10 5 cm/s, v { t 1] = 4.6 • 10 5 cm/s, v[ 2) = 3.3 • 10 5 cm/s, 

are in good agreement with the experimentally measured values vi = 6.5 • 10 5 cm/s, v t = 4.2 ■ 10 5 cm/s, v\ = 
2.6- 10 5 cm/s @.) 

Though being crude, this simple rigid-ion model gives possibility to analyze, at least qualitatively, the main proper- 
ties of the folded phonon modes in the low-temperature phase of a'-NaV205. In particular, to follow the temperature 
shifts of the normal mode frequencies at temperatures near T c , we have performed calculations in the framework of 
Fmm2 LT-structure, using the atomic coordinates from the structural work and taking into account the charge 
redistribution. Three cases have been considered, namely, (i) lattice deformation ignored, redistributed charges of 
vanadium ions in the modulated ladders, y 4 - 5±s with S — 0.05; (ii) lattice deformation ten times smaller than mea- 
sured by X-ray diffraction at 15 K in Ref. M, 8 = 0.05; (hi) the same as (ii) but 5 = 0.1. These calculations involved 
the numerical diagonalization of the dynamical matrix with dimensions 192 x 192 that was constructed using the same 
non-Coulomb force constants as in the paraphase. The main results of calculations are as follows. 

1) The soft acoustic mode at the (0, 0, ^) point of the BZ of the paraphase || becomes unstable under redistribution 
of the vanadium charge density alone, not taking into account the lattice distortion. 

2) Both lattice distortion and charge disproportion cause normal mode frequency shifts and splitting of quadruplets 
(into 2Ai + 2A2 or 2J5i + 2E>2 doublets), due to induced changes of the dynamical matrix. For different modes, 
displacements of atoms towards their new equilibrium positions and charge redistribution have effects upon frequencies 
either in the same or in the opposite direction. The shifts of "hard" singlet modes and splitting of quadruplets amount 
to 0.01-3% of their frequency. 

IV. DISCUSSION 

A. Splitting of the folded phonon modes 

First of all, we put the question whether it is possible to explain the number and polarization properties of the 
additional low-temperature modes of a'-NaX^Os supposing all of them to be vibrational folded modes of a dimerized 
lattice. Fourfold degeneracy of modes folded from the BZ Q-point at the phase transition determines some important 
peculiarities of the low-temperature vibrational spectrum of a'-NaV20s. We shall discuss them assuming, hrst, the 
Fmm2 group for the low-temperature structure |^,Q|5|. In this case, (^-quadruplets should manifest themselves as 



2Bi + 2B2 pairs of doublets having close frequencies in E || a and E || b FIR spectra. For the two narrow additional 
low-temperature modes, namely, for the lines at about 101.5 and 127 cm -1 we really observe a gradual splitting into 
such doublets upon lowering the temperature (see Figs. ||[Q,0). We suppose that the observed additional FIR lines 
with coincident or nearly coincident frequencies in E || a and E || b spectra correspond to unresolved 2B\ + 2B2 
doublets originating from Q2 quadruplets. We place them into the lower right part of Table |yL together with Raman 
data from the literature. Note, that the frequency differences between 199.0 cm -1 (E || a) and 199.5 cm -1 (E || b) 
lines and between 959.7 cm -1 (E || a) and 959.2 cm -1 (E || b) lines are beyond the error limit of our high resolution 
measurements. 

Qi-quadruplets should split into 2Ai + 2A 2 pairs of doublets. Ai modes are observed in E || c polarization, while 
A2 modes are silent in the IR spectra (but allowed in the (ab) Raman spectra). All the E || c modes are rather broad 
and we never observed a reliably resolved doublet structure. However, the shape of the lowest frequency mode near 
70 cm -1 (E || c) suggests a doublet structure (see Fig. 11). We have performed a fitting of the measured spectrum 
by the expressions (2)-(5) assuming a contribution from (i) one oscillator and (ii) two oscillators. Fig. 11 clearly 
demonstrates a much better result in the case of two oscillators. A doublet structure of the corresponding Raman 
Ai-line has been reported in Ref. |p9| . A\ and A2 Raman modes with coincident or close frequencies are listed in the 
lower left part of Table |VL So, the lower part of Table IV contains the LT modes that could originate from the BZ 



Q-point of the paraphase. In that case, 7 out of 12 Q\ and 10 out of 12 Q2 quadruplets are observed. Upper part 
of Table OV^ lists the rest of the observed additional LT modes from our measurements and from the literature. They 
could come from the Z-point, so that 10 out of 16 A\, 3 out of 6 A 2 , 11 out of 16 B lt and 3 out of 10 B 2 modes folded 
from the Z-point are observed. Raman and FIR modes that coincide in frequency within the absolute precision of 



Raman measurements (typically, ±2 cm -1 ) are put into the same line of Table IV. 

Thus, the number and polarization properties of the additional modes that appear below the temperature of the 
phase transition in a/-NaV205 do not contradict their assignment to folded modes of the distorted Fmm2 structure. 

Now we compare in more detail the experimental results to the results of lattice dynamics calculations for the 
frequency range displayed in Fig. 10. Table [V] lists the calculated frequencies of the T-point TO and LO modes and Z- 
and Q-point modes relevant to the zone folding in the case of Fmm2 LT-structure. All the experimentally observed 
modes in this frequency range can be reasonably ascribed to the modes of the Fmm2 structure (see columns 3 and 4 of 
Table 0). The modes 948 (Ai) and 948 cm -1 (A2) with coincident frequencies in Ai and A2 Raman spectra and two 
clearly different infrared modes 959.7 {B\) and 959.2 cm -1 (B2) in E || a and E || b polarizations should correspond 
to unresolved 2Ai and 2A 2 and, respectively, 2Bi and 2B 2 doublets folded from the Q-point. According to our lattice 
dynamics calculations (for the case (hi) of the Section IIIC), Q\ (969 cm -1 ) and Q2 (973 cm -1 ) quadruplets (see 
Fig. 10) split into 2A X (967.8-970.3 cm" 1 ) and 2A 2 (969.3-969.4 cm" 1 ) and, respectively, 2B X (974.8-977.0 cm" 1 ) 
and 2B2 (973.5-973.6 cm -1 ) doublets. These splittings are within the limits of the linewidths. 

If the space group for the LT-phase of a'-NaV205 were a subgroup of Fmm2 (e.g., Ccc2 |lj],[l8|] ) , 48 additional 
doublets 16 x 2Ai + 6 x 2A2 + 16 x 2B\ + 10 x 2B2 from the A point would appear independently in the spectra 
of different symmetries and also one would expect additional 96 singlets 24(Ai + A2) and 2A(B\ + B2) with close 
frequencies in A\ and A2 and B\ and B2 spectra, originating from the folded S- and i?-doublets. Domains due to 



oblique charge ordering pattern possible for a zigzag charge order in every ladder in the space group Ccc2 [117 32 
would not affect the direction of the orthorhombic a- and b-axes and, thus, polarization properties of vibrational 
modes. Experimental data on the new LT-modes do not contradict both Fmm2 and Ccc2 symmetries. However, 
while in the first case one has to admit that 27 new frequencies remained unobserved, in the latter one this quantity 
grows up to 123. In both cases, close or coincident frequencies in A\ and A2 and B\ and B2 spectra are naturally 
explained by their origin from the folded BZ Q-point (Fmm2) or Q, A, S, and i?-points (Ccc2). A possible monoclinic 
distortion (space groups C2 pll or C2/c pi) would violate the E || a or E || b FIR selection rules, as well as Raman 
selection rules. Probably, the growing contribution of the Ai g modes 530 and 304 cm -1 to the (ab) Raman spectrum 
for the temperatures below ~ 100 K observed in Ref. p8f is just because of such a distortion. 

At the first glance, LT modes fit the C2/c centrosymmetric group proposed in Ref. jl9|. In particular, in the case 
of C2/c group (^-quadruplets would split into 2B U (E || a, E || b) + 2B g (ac,bc) doublets. Such a scenario seems 
to describe better the experimental data (see the lower right part of Table Evfi. However, the frequency differences 
between the observed FIR and Raman modes are within the precision of Raman measurements. 

To summarize, it is not possible to chose between the proposed space symmetry groups for the low-temperature 
crystal structure of a'-NaV20s on purely spectroscopic grounds. The number and polarization properties of all the 
observed new LT modes in a'-NaV20s can be reasonably described within the conception of folded vibrational modes 
of the dimerized Fmm2 crystal structure. Observed FIR doublets and close frequencies in different polarizations 
are naturally explained by their origin from the Q-point modes folded into the zone center. Prior to analyse the 
temperature dependences of frequencies and widths of new modes and to discuss alternative interpretations, we note 



that published FIR [E0 23 2qJ5U32fl and Raman p,El,E7y2a,iqBaE9 data agree with each other, within the precision 



of measurements (typically, ±0.5 cm for FIR and ±1.5 cm for Raman measurements). 



B. Magnetic bound states? 

As we mentioned in Section IIIA, three of the observed new low-temperature modes, namely, modes near 70, 107 
and 133 cm -1 in E || c polarization, have different frequencies in different samples (see Table IV). These modes and 



the mode near 125 cm 1 (E || c) demonstrate unusually large broadening, shift, and loss of intensity upon heating (see 



Figs. 6, 8). All these modes were observed also in Raman scattering (see Table IV), exhibiting the same properties 
as their infrared counterparts. 

In Refs. [pitp9[ , Raman modes near 65 and 132 cm -1 have been assigned to one- and two-magnon scattering 
processes, respectively. The frequency 64-70 cm -1 of the lowest- frequency sharp feature in FIR and Raman spectra 
below T c corresponds well to the size of the gap Ao = 8.13 meV= 65.6 cm -1 at k = obtained from submillimeter 
ESR measurements |37J • The transitions between the singlet ground state and the excited triplet state are forbidden 
(AS* = 1). The Dzyaloshinsky-Moria (DM) interaction with DM vector along the c-axis breaks the selection rule for 
ESR (and FIR) E || c transitions in a'-NaV 2 5 f^Q. Recently, it has been shown theoretically @ that for a'- 
NaV2C>5 the DM interaction results in finite Raman cross-section for both one- and two-magnon scattering processes. 
A single-magnon Raman line of this type should show no splitting in an external magnetic field parallel to the c-axis 
but should split into two branches for a field perpendicular to the c-axis. The authors of Ref. p3] communicated 
that the 67 cm -1 Raman line neither shifts nor splits nor even broadens in a magnetic field up to 7 T. They have 
attributed this line and the 107 and 134 cm -1 Raman modes to magnetic singlet bound states, on the basis of (i) 
the temperature dependence of the intensity, frequency shift, and lincwidth, and (ii) the frequency dependence on a 
sample quality. All the three modes rapidly loose their intensity, soften and broaden on approaching T c from below 
pTLESfl . It is worth mentioning that very recently two more Raman modes with similar properties have been reported 
(namely, 86 and 126 cm -1 modes in (ab) geometry) |29| . To our opinion, the interpretation of the 67, 107, and 
134 cm -1 modes as magnetic bound states does not follow unambiguously from the above mentioned arguments. 

One of the main arguments in favor of considering 67, 107, and 134 cm -1 Raman modes as magnetic bound states 
was a slow increase of their intensity upon cooling. For comparison, the intensities of the folded phonon Raman modes 
at 202, 246, and 948 cm -1 rose sharply below T c and then saturated. However, our measurements show that many 
of the transition- induced modes demonstrate a slow growth of intensity upon cooling below T c . Among them, the 
considered mode near 107 cm -1 (E || c) and the folded phonon mode 111.7 cm" 1 (E || b) behave in the same way 
(see Fig. 8). 

Different temperature dependences of the type (T c — T) 2n " (here j3 is the critical index of the order parameter; 
n = 1,2,3) for the intensities, frequencies and widths of the folded modes in the vicinity of a structural phase 
transition have been reported and analysed by Petzelt and Dvorak EG], taking into account anharmonic couplings 
between modes. In a'-NaV2C>5, as we show below, a mere charge disproportion may lead to the (T c — T) 2n ^ 1 n > 1, 
intensity dependence for a folded phonon mode below T c , even in the harmonic approximation. 

Intensity of a FIR spectral line is proportional to a square of the electric dipole moment of a primitive cell induced 
by the respective vibrational mode. The dipole moment may be represented as a sum of products of effective atomic 
charges and corresponding displacements of Bravais sublattices. A slow intensity dependence may be expected when 
the main contribution to the dipole moment comes from the charge disproportion terms that involve displacements 
renormalized by the order parameter. Let us consider, for example, the z-component of the electric dipole moment 
due to displacements of vanadium ions that transform according to the A\ representation of the factor group of the 
LT phase. The primitive cell of the Fmm2 structure contains eight V ions that form two rectangular plaquctts 
(fragments of nonmodulated ladders) in the neighboring afe-layers, four V 4+ , and four V 5+ ions on the rungs of 
the neighboring (left and right) modulated ladders (see Fig. pi). Normal modes that correspond to Z\, Z^, and Q\ 
representations involve antiphase translations of V 45+ plaquetts along the c-axis and antiphase displacements of V 5+ 
and V 4+ ions along the same axis. At T c , all the displacements have equal amplitudes, vanadium ions have equal 
charges, and the corresponding electric dipole moment is zero. At T < T c , nonzero dipole moment appears due to 
different amplitudes of the plaquette translations and different amplitudes of V 45_(5 and V 45+(5 ions. Our calculations 
confirmed that, for some low frequency modes involving large amplitudes of vanadium ions vibrations, the effective 
electric dipole moment is determined mainly by the terms containing the charge disproportion 8 and the differences 
between ion amplitudes at the temperatures T and T c . These terms provide the (T c — T) 4/3 dependence of the FIR 
absorption. The same consideration is applicable to the Raman spectra where intensities are determined by electric 
polarizabilities of lattice modes. 

Thus, a slow growth of the mode intensity below T c is not a unique property of the spin-gap related excitations 
only. Futhermore, if we compare the intensities of the FIR new modes measured in this work with the integrated 
intensity of the spin-gap mode observed in optical ESR measurements |37[| , it is evident that only the 133 cm -1 mode 
follows the same temperature dependence as the gap-mode (see Fig. 8). 

Now in Fig. 12 we compare the normalized magnetic gap energy, as observed in neutron scattering experiments pOj, 



and the normalized frequencies of the 70 and 106.9 cm -1 FIR modes. The latter change much more slowly on 
approaching T c than the magnetic gap. We could not measure the frequency of the weak 133 cm -1 FIR mode below 
24 K but the corresponding Raman mode practically does not shift within the interval of T/T c from 0.2 to 0.88, as one 
can see from the spectra of Rcf. [g9|. Thus, the temperature dependence of the peak positions for the lines near 70, 
107, and 124.7 cm -1 differ significantly from the dependence for the singlet-triplet gap, A (T), while the frequency of 
the Raman peak near 30 cm -1 corresponding to a magnetic bound state in CuGeC>3 followed the same temperature 
dependence as A (T) Q. 

To summarize, the comparison between the temperature-dependent intensities of the transition-induced FIR modes 
and of the spin-gap mode show that the mode 133 cm -1 is the only one that could be considered as a magnetic bound 
state. However, the temperature dependence of its frequency contradicts this hypotesis. Thus, the assignment of the 
70, 107, and 133 cm -1 modes to magnetic bound states seems to be wrong. 

As an alternative interpretation we would propose folded phonon modes that strongly interact with charge and 
spin excitations. The interaction renormalizes the frequencies of phonons, causes mode broadening and softening 
near T c . A sodium deficiency strongly influences magnetic (and charge) states which leads to a strong change of 
interaction parameters. This circumstance could explain the dependence of 70, 106, and 133 cm -1 mode frequencies 
on a sample quality. Our earlier observation on the ir-dependent phonon frequencies in Nai_j;V205 la] supports such 
an assumption. We observed a giant frequency shift of 30 cm -1 for one of the (aa) Raman modes. Namely, the 
447 cm -1 mode [x = 0) moved to 477 cm -1 for x = 0.15. Simultaneously, the maximum of a broad-band scattering 
moved in the opposite direction, from 632 cm" 1 for x — to 562 cm" 1 for x = 0.15. It was naturally to explain this 
change of the frequency difference between two Raman bands by a change of interaction between the corresponding 
excitations ||. Though the exact nature of the broad band in the (aa) Raman geometry is not yet clear, it must be 
connected with low-frequency electronic excitations on the ladder rungs |28| , p5| . 

The interaction between spin, charge, and lattice degrees of freedom is an intrinsic property of a'-NaV20s . In sodium 
vanadate one d-electron belongs to two vanadium atoms which gives an additional electronic degree of freedom. The 
electron hopping amplitude along a ladder rung is significantly higher than the amplitudes of hopping along a ladder 
or between different ladders []2|,|s| . Zigzag charge ordering at low temperatures can be described as an antiferroelectric 
ordering of electric dipoles on the rungs of vanadium ladders pa]. As each electron carries the spin (S=l/2), charge 
and spin degrees of freedom are closely connected. To describe this situation, Mostovoy and Khomskii proposed a 
spin-isospin model and used it to show that the antiferroelectric charge ordering in the system of vanadium ladders 
opens a spin gap J42j,|35j. They argued that in a'-NaV2C>5 the quasi-one-dimensional spin system is strongly coupled to 
a low-energy antiferroelectric mode (isospin mode) of the excitonic type and the observed charge ordering corresponds 
to the softening of this mode. The softening of the isospin excitation occurs at the wave vector of the superlattice 
structure appearing in the ordered phase p5j, that is, at the Q-point {\\\)- Low-energy vibrational excitations of 
the corresponding symmetry should strongly interact with the isospin mode. To understand a possible mechanism 
of such an interaction we refer to Fig. 13 showing atomic displacements for vanadium atoms in the lowest-frequency 
A\ mode (at 136 cm -1 , according to our calculations) folded from the BZ Q-point into the zone center. Vibrations 
are mainly along the c-axis, vanadium atoms at the two ends of a rung move in antiphase, that is when the left 
atom approaches the apical oxygen, the right one moves away from it. Such a motion stimulates vanadium charge 
redistribution (V 45±<5 ) jl3|] that corresponds just to a static picture of charge ordering. 

By inspecting the experimental results summarized in Table | we notice that the lowest frequency ^li-symmetry 
modes (E || c) are about an order of magnitude broader than the lowest-frequency B\ (E || a) and B2 (E || b) modes 
even in the limit T — » 0. They have the ratio j/u> ~ 2 ~ 6 • 10" 2 in comparison with 1-^3- 10" 3 for B\ and B2 
modes. Such a broadening could result from a decay of the low-frequency vibrations into isospin and spin excitations, 
due to the interaction with them. High-resolution neutron inelastic scattering experiments on the spin excitations of 
a'-NaV2C>5 revealed the two branches associated with distinct energy gaps QJ45|. At the zone boundary (||0) (as 
well as at the zone center) these gaps constitute 68.6 cm" 1 Q (67 cm" 1 according to Ref. Q) and 91.1 cm" 1 |44|] 
(87.8 cm" 1 J45[). Keeping in mind the energy and momentum conservation laws, the following decay channels are 
possible: direct interaction with the zone center magnon or a decay into two zone boundary isospin excitations (for 
the 70 cm" 1 mode), a decay into a zone boundary magnon (66 or 90 cm -1 ) and the isospin excitation (for the 106 or 
125 cm -1 modes). Such a process corresponds to a hopping of an electron along the ladder rung and a simultaneous 
spin-flip. In that case, the zone boundary isospin gap would constitute of about 35 cm" 1 , in reasonable agreement 
with the theoretical estimates [tS5|. 



C. In search for a soft mode 

The lattice dynamics calculations for the high-temperature phase of a'-NaV20s predicted a soft mode behavior of 
the transverse acoustic mode at the BZ boundary with k = (00|). It was found to become unstable under charge 
redistribution below T c . So, one could expect a doubling of the unit cell along the c-axis, as a precursor of the 
subsequent magnetic ordering with a dimerization along the a-, b-, and c-directions. We have to mention that two 
subsequent phase transitions within the temperature interval of 1 K have been communicated in Ref. [|46| . 

The discussed soft mode corresponds to the irreducible representation Z5 of the BZ (00 5) point and it could be 
observed as a folded FIR mode in E || a polarization. We have thoroughly studied low-temperature a-polarized FIR 
transmittance spectra of 1mm x 3mm x 111/xm sample of a'-NaX^Os in the spectral range 20-138 cm' 1 but failed 
to find a soft mode. Possibly, it has a frequency lower than 20 cm" 1 at 6 K or its oscillator strength is too small. 
Another possibility would be that this soft mode is simply a property of our rather crude model of lattice dynamics. 

The E || c polarized 70 and 107 cm -1 modes demonstrate the strongest observed softening toward T c . But this 
softening is only of about 3 % and is, probably, connected with the interaction of lattice modes with magnetic and 
charge excitations, as discussed in the previous section. It looks like the true soft mode of the phase transition in 
a'-NaV205 is of electronic nature (low-energy antiferroelectric (isospin) mode [p5fl). 

V. CONCLUSIONS 

In summary, we report the temperature-dependent polarized high-resolution FIR transmission spectra of a 1 - 
NaV2C>5. The comparison with the results of neutron scattering and optical ESR measurements and the theoret- 
ical analysis of the properties of vibrational modes in a charge ordered state lead us to the conclusion that all the 
transition-induced new modes of a'-NaX^Os can be completely described in terms of folded modes of a dimerized 
lattice. 

Based on the present results, the widely accepted interpretation of several lowest-frequency new modes as magnetic 
bound states seems highly unlikely. As an alternative interpretation of 70, 106, 125, and 133 cm -1 spectral lines 
we suggest low-frequency lattice vibrations of an appropriate symmetry that strongly interact with charge and spin 
excitations at the BZ boundary of the high-temperature phase and fold into the zone center at the phase transition 
into a dimerized state. 

Note: When this paper already was with editors be became aware of the work J17J where E || c transmittance 
of a'-NaX^Os was measured and low-frequency transition-induced modes were reported to appear at 68, 106, and 
124 cm -1 (T — 4.2 K) and not to change under the magnetic field up to 15 T. 
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TABLE I. Parameters at T — 9 K of the FIR absorption lines activated by the phase transition in «'-NaV2 05. The 
frequencies of E || c modes in the sample 2 are given in brackets. Here uij is the mode frequency, 7j is the damping constant, 
fj is the oscillator strength. 





E||a 






E || b 






E||c 




ujj, cm 1 


7j> cm 1 


io 4 ■ u 


uij, cm 1 


7j, cm -1 


10 4 • u 


ujj, cm 


Ij, cm 4 


10 4 • fj 


91.2 


0.17 


4.7 


101.4 


0.20 


5.2 


70.0 (67.5) 


4.5 


33 


101.4 


0.24 


8.0 


111.7 


0.24 


1.2 


106.9 (106.2) 


2.3 


10 


101.7 


0.16 


3.6 


126.8 


0.31 


3.4 


124.7 (124.7) 


1.0 


1.6 


125.6 


0.24 


0.4 


127.5 


0.33 


9.0 


133.1 (131.6) 


2.5 


0.5 


126.8 


0.33 


13.0 


199.5 


1.3 


2.5 


140.0 


0.7 


1.3 


127.6 


0.30 


4.8 


234.2 


2.0 


2.6 


148.2 


0.6 


0.7 


145.0 


0.46 


3.2 


324.6 


3.9 


5.2 


256.5 


2.6 


0.05 


145.7 


0.43 


7.6 


410.3 


2.0 


1.5 


327.4 


6.0 


0.17 


147.9 


0.5 


26 


959.2 


1.3 


0.06 








157.2 


0.3 


0.8 














199.0 


1.8 


46.2 














410.4 


4.3 


26 














451.1 


2.9 


3.7 














959.7 


1.5 


0.4 















TABLE II. Wyckoff positions of atoms in the elementary cells of a'-NaV20s above and below T c . Figures before chemical 
symbols denote numbers of atoms of the same kind. 





Atom 


T > T c , Pmmn g 
Wyckoff position 


Site 


symmetry 




Atom 


T < T c , Fmm2 § 
Wyckoff position 


Site 


symmetry 


1 


V 


4/ 






cr 


2 

4 




VI 
V2 




16e 
8d 






C x 

cr 


1 


Na 


26 






C*2i, 


4 
2 




Nal 
Na2 




4a 
86 






Civ 

c 2 


1 


Ol 


2a 






C?v 


2 
2 




Oil 
012 




8c 
8d 






cr 
cr 


1 


02 


4/ 






cr 


2 
4 




021 
022 




16e 
8d 






Ci 

cr 


1 


03 


4/ 






cr 


2 

4 




031 
032 




16e 
8d 






d 
cr 



TABLE III. Number of atoms in a primitive cell n, number of BZ center optical vibrational modes N; points in the BZ of 
the paraphase that fold into the zone center and corresponding optical modes in the cases of different proposed space groups 
for the distorted structure, m is the number of additional frequencies just below T c . 



Space group 


n 


N 


BZ points 


Optical modes 


m 


Fmm2 [|0JI|] 


64 


189 


r, z, q 


45r + 48Z + 24 x 4Q 


72 


Ccc2 F7 
C2/c |1£ 
andC2 


1 


128 


381 


F, Z, Q, A, S, R 


45r + 48Z + 24 x 4Q 
+48 x 2A + 24 x 25 
+24 x 2R 


168 
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TABLE IV. Additional low-temperature IR and Raman modes of «'-NaV20s. 
used in this work or in Ref. |2q|. 



Figures in brackets refer to the second sample 





Ai 


A 2 


Bx 


B 2 


IR 


Raman 


Raman 


IR 


Raman 


IR 


Raman 


E||c 


(aa), (bb), (cc) 


(ab) 


E || a 


(ac) 


E || b 


(be) 


107 (106) 


107 (103) a 


86 c 


91.2 


116 e 


111.7 




140.0 


140 6 




125.6 




234.2 






164 Q 


159 6 , 158 c 


145.1 






354 e 


* 


202" 


422 a,d 


145.7 








256 


296 6 




147.9 
157.2 


157 J 






314 






* 


306^ 






* 


394 a 




374 s 








* 


650 a 




451.1 








* 


692 Q 




787.3 
[969*] 








70.0 (67.5) 


67 (64) a 


67 (64) a 


101.4 
101.7 


102 / 


101.4 




124.7 (124.7) 


127 6 , 123' 


127 6 , 126 c 


126.8 




126.8 




133.1 (131.6) 


134 (131) a 


134 (131) a 


127.6 
199.0 


128 7 
200 / 


127.5 
199.5 




148.2 


151 a 


149 6 , 148 c 


* 


323 e 


324.6 


323 e 


244 


246 a 


24&a ,d 


362 s 




362 s 




327.4 


325 a 


332 a,d 


410.3 


408° 


410.3 


408 e 


* 


948 a 


948 a 


* 


433 e 


* 


433 e 








717 s 


720 e 


718 s 


720 e 








959.7 




959.2 





* Strong absorption due to HT-phonons in this frequency region 

a Low-temperature Raman data from Ref. [E8[ 

b Weak feature in the low-temperature spectra of Ref. 1 28 

c Low-temperature Raman data from Ref. p9] 

d Fano resonance 

e Low-temperature Raman data from Ref. pll 

' Weak feature in the low-temperature spectra of Ref. |L9[ 

9 Low-temperature FIR data from Ref. J32J 

Nonpolarized FIR transmittance data for a powdered sample from Ref. |26 
1 Weak feature in the low-temperature (aa) spectrum of Ref. [g9[ . 
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TABLE V. Comparison between the calculated vibrational frequencies of the F, 
Pmmn structure and measured frequencies for the HT- and LT-phases of a'-NaV^Os 
for the Fmm2 LT-structure are indicated. 



Q, and Z-modes for the high-temperature 
Corresponding irreducible representations 



Pmmn 




calculated 




measured 


Fmm2 


r 


B 3 u 


(TO-LO) 


955-957 


939 (E || a) | 


Si 




Bin 


(TO-LO) 


961-1036 


955-1017 (E || c) § 


Ai 




B 2 g 




962 


951 (ac) § 


Si 




A g 




964 


970 (cc, aa) § 


A 1 


Q 


Qi 




969 


948 (aa, bb)-948 (ab) §| 


2Ai + 2A 2 




Q 2 




973 


959.7 (E || a)-959.2 (E || b) 


2Bi + 2B 2 


z 


z 5 




985 


969 a @ 


Bi 




z 6 




1003 


* 


Si 




Zi 




1018 


* 


A 1 




z 5 




1024 


* 


A 1 



Strong absorption due to HT-phonons in this frequency region 
Observed in unpolarized transmittance spectra of a powder sample 
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FIGURE CAPTIONS 

FIG. 1. The structure of the low-temperature phase of a'-NaA^Os. (a) Grey pyramids incorporate V 4 ' 5+ ions. Black (white) 
pyramids are bigger (smaller) in size and are occupied by the V 4+ (V 5+ ) ions. Balls represent Na atoms, (b) The ab-projection 
of the one V-0 layer. Black (white) circles stand for V 4+ (V 5+ ), while the grey ones represent V 4 ' 5+ . Arrows indicate the 
direction of displacements at the phase transition into the LT-phase. 

FIG. 2. Transmittance spectra of a'-NaV 2 5 at T « 40 K> T c (dashed line) and T = 8 K< T c (solid line) in E || a, E || b, 
and E || c polarizations of the incident light with k || c (a,b) and k || a (c). The sample thicknesses were 14 ^m (a,b) and 
150 fim (c). Arrows show new low-temperature modes below T c . 

FIG. 3. Absorbance difference spectra in E || a polarization for different temperatures showing the growth of new spectral 
lines below T c . Spectra are vertically shifted for the sake of clarity. 

FIG. 4. Absorbance difference spectra in E || b polarization: transition-induced modes with frequencies below 230 cm" 1 . 

FIG. 5. Absorbance difference spectra in E || c polarization. An asterisk marks weak features due to a leakage of E || a 
polarization. A weak intrinsic c-polarized mode near 133 cm -1 is shown by an arrow. 

FIG. 6. Temperature dependences of the normalized shifts and linewidths of transition-induced modes. The insets show 
temperature dependences for the two modes from the high-temperature phase of a'-NaV2 0s. 

FIG. 7. Peak position of selected transition-induced modes at different temperatures. 

FIG. 8. Temperature dependence of the normalized oscillator strengths /t//sk for several transition-induced modes of 
a'-NaV^Os. Lines are guides for the eye. Intensities of the X-ray reflections from Ref. [|10| are shown by grey symbols. 
Asterisks represent the normalized integrated intensity of the spin-gap mode observed in optical ESR, as estimated from the 
data of Ref. Q. 

FIG. 9. Brillouin zone of the Pmmn HT-structure of a'-NaV20s. 

FIG. 10. Calculated phonon dispersion for the high-temperature Pmmn structure of a'-NaV^Os. 

FIG. 11. Lowest-frequency transition-induced mode near 70 cm" 1 in E || c transmittance spectrum of the 150 /im thick sam- 
ple, corrected for the additional absorption of the epoxy layer (open circles). The dashed line was calculated using Eqs. (2)-(5) 
and one-oscillator approximation {uj — 70.0 cm" 1 , 7 = 4.5 cm" 1 , / = 3.3 ■ 10" 3 ). The solid line represents the fit result for the 
two oscillators (wi = 68.75 cm" 1 , u 2 = 70.9 cm" 1 ; 71 = 2.2 cm" 1 , 72 = 3.2 cm" 1 ; /1 = 1.2 • 10~ 3 , f 2 = 1.85 ■ 10" 3 ). 

FIG. 12. Comparison between temperature dependences of normalized energies (frequencies) fUVr/Ttuis k of the transi- 
tion-induced modes near 70 and 107 cm -1 and normalized magnetic gap At/A7 k form inelastic neutron scattering measure- 
ments llol (black circles). 
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FIG. 13. Picture of atomic displacements of vanadium atoms for the lowest-frequency A\ mode (Fmm2) folded from the 
Q-point of (Pmrnn) high-temperature structure, (a) one ab-layer. Only vanadium atoms are shown. Displacements along the 
positive (negative) direction of the c-axis are indicated by points (crosses), (b) Vanadium atoms 1 and 2 at the opposite ends 
of a ladder rung, together with oxygen at the middle of a rung and two apical oxygens. When V2 approaches apical oxygen, 
VI goes moves away from it, which stimulates an asymmetric charge distribution along the rung. 
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